(12) DEMANDE INTERNATIONALE PUBLIEE EN VERTU DU TRAITE DE COOPERATION 

EN MATIERE DE BREVETS (PCT) 



(19) Organisation Mondiale de la Propriete 
Intellectuelle 

Bureau international 

(43) Date de la publication internationale 
12 decembre 2002 (12.12.2002) 



PCT 



llllllllllllllllllllllllllllllllllllllllllllllll 

(10) Numero de publication internationale 

wo 02/100111 Al 



(51) Classification internationale des brevets 7 : H04N 7/26, 
G06T 9/00, 17/20 



Rennes (FR). LAURENT-CHATENET, Nathalie 

[FR/FR]; 3, rue des Fraiches, F-35690 Vignoc (FR). 



(21) Numero de la demande internationale : 

PCT/FR02/01937 



(74) Mandataire : GUENE, Patrick; Le Nobel, 2, allee An- 
toine Becquerel, BP 90333, F-35703 Rennes Cedex 7 (FR). 



(22) Date de depot international : 6 juin 2002 (06.06.2002) (81) Etats designes (national) : AE, AG, AL, AM, AT, ALT, AZ, 

BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, CZ, 
francais DE, DK, DM, DZ, EC, EE, ES, FT, GB, GD, GE, GH, GM, 

HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, LK, 
francais LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, MX, 

MZ, NO, NZ, OM, PH, PL, PT, RO, RU, SD, SE, SG, ST, 
SK, SL, TJ, TM, TN, TR, TT, TZ, UA, UG, US, UZ, VN, 
YU, ZA, ZM, ZW. 



(25) Langue de depot : 

(26) Langue de publication : 



(30) Donnees relatives a la priorite : 

01/07406 6 juin 2001 (06.06.2001) FR 



(71) Deposant (pour tons les Etats designes sauf US) : (84) Etats designes (regional) : brevet ARIPO (GH, GM, KE, 



FRANCE TELECOM 

F-75015 Paris (FR). 



[FR/FR]; 6, Place d'Alleray, 



(72) Inventeurs; et 

(75) Inventeurs/Deposants (pour US seulement) : DE- 
MARET, Laurent [FR/FR]; 32, rue Francis Toly, F-35000 



LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), brevet 
eurasien (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), brevet 
europeen (AT, BE, CH, CY, DE, DK, ES, FT, FR, GB, GR, 
IE, IT, LU, MC, NL, PT, SE, TR), brevet OAPI (BE, BJ, 
CF, CG, CI, CM, GA, GN, GQ, GW, ML, MR, NE, SN, 
TD, TG). 

[Suite sur la page suivante ] 



(54) Title: METHODS AND DEVICES FOR ENCODING AND DECODING IMAGES USING NESTED MESHES, 
PROGRAMME, SIGNAL AND CORRESPONDING USES 

(54) Titre : PROCEDES ET DISPOSITIFS DE CODAGE ET DE DECODAGE D' IMAGES METTANT EN OEUVRE DES 
MAILLAGES EMBOITES, PROGRAMME, SIGNAL ET APPLICATIONS CORRESPONDANTES 




(57) Abstract: The invention concerns a method for encoding at least a 
source image using a hierarchical mesh defining at least two nested spaces 
corresponding each to a decomposition level of said meshing, which 
consists, at least at a decomposition level n (except the first decomposition 
level), in delivering only image coefficients expressed in a base of functions 
defined in a space orthogonal to the space corresponding to the preceding 
decomposition level n-1, said functions being selected so that said image 
coefficients enable to optimise for said decomposition level n the data 
already transmitted for the preceding decomposition level n-1, so as to 
produce a reconstructed image, representing said source image, with 
optimised restoration quality for said decomposition level n. 

(57) Abrege : L' invention concerne un procede de codage d'au moins 
unc image source mcttant en ocuvrc un maillagc hicrarchiquc definissant 
au moins deux espaces emboites, correspondant chacun a un niveau de 
decomposition dudit maillage, selon lequel, a au moins un niveau n de 
decomposition (a 1' exception du premier niveau de decomposition), on 
delivre uniquement des coefficients d' image exprimes dans une base de 
fonctions definic dans un cspacc orthogonal a l'cspacc correspondant au 
niveau precedent n-1 de decomposition, lesdites fonctions etant choisies 
de facon que lesdits coefficients d' image permettent d'optimiser pour 
ledit niveau n de decomposition les informations dejo transmises pour le 
niveau precedent n-1 de decomposition, de facon a produire une image 
reconstruite, representative de ladite image source, avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 
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Procedes et dispositifs de codage et de decodage damages mettant en 
ceuvre des maillages emboites, programme, signal et application 
correspondantes. 

Le domaine de 1' invention est celui du codage d' images, en vue 
5 notamment de leur transmission ou de leur stockage. Plus precisement, l'invention 
conceme une amelioration des techniques de codages hierarchiques, mettant en 
oeuvre une hierarchie de maillages emboites. 

L'approche des maillages emboites a deja fait l'objet de nombreuses 
etudes. Elle est par exemple presentee dans [7] (par souci de lisibilite, les 
10 differents documents cites sont regroupes en annexe 3, a la fin de la presente 
description), dans le cas du codage d' images video. 

Un maillage est classiquement defini par un ensemble de sommets et de 
faces orientees (figure 1) definissant une topologie. De tels maillages sont par 
exemple utilises en graphisme sur ordinateur, pour modeliser des objets en trois 
15 dimensions avec une complexity geometrique limitee. 

L' approximation d'un maillage M consiste a trouver un maillage M' dont 
la complexity geometrique est moindre que celle du maillage M, et qui approche 
au mieux la geometrie de M. 

Avantageusement, le maillage M' est constitue d'une succession de 
20 maillage emboites, correspondant chacun a un niveau de detail, ou niveau 
hierarchique, de fagon a permettre une reconstruction progressive des images et 
un codage simplifie. 

A chaque niveau hierarchique, les valeurs nodales du maillage sont 
optimisees pour minimiser Terreur quadratique de reconstruction. Ces valeurs 
25 nodales sont ensuite quantifies et codees. Une telle methode permet d'atteindre 
des taux de compression performants et de limiter les degradations visuelles qui 
correspondent plus ici a des effets de lissage, moins desagreables pour Pceil 
humain. Ceci est lie aux bonnes proprietes de continuity des surfaces reconstruites 
grace a V approximation par les maillages. 
30 De plus ce schema s'avere adapte a la video. En effet, les maillages 



WO 02/100111 



2 



PCT/FR02/01937 



triangulares s'averent plus souples et efficaces pour 1' estimation du mouvement. 

Les inventeurs ont cependant constate que cette technique de compression 
souffre d'un defaut de sous-optimalite. En effet, les bases utilisees pour chaque 
niveau presentent des redondances. 
5 En outre, la technique des maillages emboites permet de reconstruire 

progressivement des images, tout d'abord d'un niveau grossier, puis 
progressivement raffine (transmission "scalable" des images). Selon la technique 
classique, cette approche n'est cependant pas optimisee : en effet, la qualite 
obtenue a chaque niveau (a 1' exception du dernier) n'est pas optimale. 

10 On connait par ailleurs des techniques de compression d'images basee sur 

F exploitation des correlations entre sous-bandes issues d'un filtrage mettant en 
oeuvre des ondelettes. Ces techniques autorisent des taux de compression tres 
performants [4] [5]. 

Cependant pour des applications pour lesquelles on ne dispose que d'un 

15 debit tres bas, on observe pour des methodes une degradation visuelle non 
negligeable en particulier sous la forme d'effets oscillatoires le long des contours. 

Par ailleurs les bases d 5 ondelettes utilisees sont construites comme produit 
tensoriel d' ondelettes mono-dimensionnelles. Ceci induit une limitation de la 
capacite a representer certaines structures en privilegiant certaines directions. Ces 

20 defauts propres aux methodes de codage d'images par sous-bandes ont conduit h 
privilegier une representation de F image basee sur des maillages triangulaires 
emboites. 

Dans [1], E.Quak suggere d f utiliser simultanement les deux techniques, en 
associant a chaque niveau du maillage une base d' ondelettes complementaires en 
25 donnant des conditions sur les aretes. II construit ainsi une base explicite de 
preondelettes sur un maillage triangulaire. Cette technique est destinee a la 
representation et a la compression de modeles numeriques de terrains en 3D. 
On presente plus en detail en Annexe 1 la maniere dont Quak precede. 
L' Annexe 2 rappelle quant a elle le principe general et les grandes lignes 
30 du fonctionnement de la methode de codage basee sur une hierarchie de maillages 
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Bien qu'efficace, la technique connue decrite dans cette Annexe 2 presente 
certaines limites. 

Notamment, un inconvenient de cette methode, telle que proposee dans 
5 [7], reside dans la non-orthogonalite des fonctions (pi (p) avec l'espace vectoriel Vp^ 
engendre par les fonctions tp i (p ~ 1) . Cela implique une concentration d'energie 
moins efficace que pour les transformees utilisant des bases de transformation 
orthogonales. 

De plus, cela pose un probleme de cohabitation de differents niveaux de 
10 resolution. Ainsi lorsqu'un sommet se trouve a la frontiere entre une zone raffinee 
et une zone non raffinee on ne peut choisir de valeur optimale pour ce sommet 
pour les deux resolutions. En effet, pour un tel sommet choisir une des valeurs 
d'un des niveaux d' optimisation ne permettra d' avoir la reconstruction optimale 
que sur la region correspondant a ce niveau de raffinement. 
15 Par ailleurs cette sous-optimalite de la representation constitue aussi un 

inconvenient dans le cadre d'un schema de codage scalable. En effet, elle ne 
permet pas de fournir la qualite de reconstruction optimale pour des debits 
intermediates. 

L'invention a notamment pour objectif de pallier ces inconvenients de 
20 l'etatdel'art. 

Plus precisement, un objectif de l'invention est de fournir une technique de 
codage d' image scalable et une technique de decodage correspondant permettant 
d'obtenir une qualite de reconstruction optimale a chaque niveau de 
reconstruction. 

25 Un autre objectif de l'invention est de fournir de telles techniques de 

codage et de decodage qui necessitent un debit limite, pour chaque niveau de 
reconstruction. 

L'invention a egalement pour objectif de fournir de telles techniques de 
codage et de decodage, qui permettent de traiter efficacement plusieurs images 
30 ayant la merae structure (meme taille et meme maillage de reference). 



WO 02/100111 



4 



PCT/FR02/01937 



Encore un autre objectif de l'invention est de fournir une structure de 
donnees et de signal qui permette d'optimiser le debit necessaire pour la 
transmission et le stockage d'images codees de la sorte. 

Ces objectif s, ainsi que d'autres qui apparaitront plus clairement par la 
5 suite, sont atteints a l'aide d'un procede de codage d'au moins une image source 
mettant en oeuvre un maillage hierarchique definissant au moins deux espaces 
emboites, correspondant chacun a un niveau n de decomposition dudit maillage. 

Selon l'invention, a au moins niveau n de decomposition (a l'exception du 
premier niveau de decomposition), on delivre uniquement des coefficients 
10 d'image exprimes dans une base de fonctions definie dans un espace orthogonal a 
l'espace correspondant au niveau precedent n-1 de decomposition, lesdites 
fonctions etant choisies de fagon que lesdits coefficients d'image permettent 
d'optimiser pour ledit niveau n de decomposition les informations deja transmises 
pour le niveau precedent n-1 de decomposition, de fagon a produire une image 
15 reconstruite, representative de ladite image source, avec une qualite de restitution 
optimisee pour ledit niveau n de decomposition. 

Avantageusement, a tous lesdits niveaux de decomposition (a l'exception 
dudit premier niveau de decomposition), on delivre uniquement des coefficients 
d'image exprimes dans une base de fonctions definie dans un espace orthogonal a 
20 l'espace correspondant au niveau precedent n-1 de decomposition. 

De fafon avantageuse, on associe audit niveau n de decomposition un 
espace mathematique de representation de Timage V n tel que V n = V n-1 ® W n . l5 
ou : 

- V n _! est l'espace associe au niveau n-1 de decomposition ; 
25 - W n _! est un espace orthogonal a V n . l5 

et on associe audit espace V n d'une part une base de fonctions affines par 
morceaux cp n , et d' autre part une base formee de la combinaison de deux bases 
orthogonales : 

- une base de fonctions affines par morceaux (p aA dudit espace V n4 ; 
30 - une base de fonctions de pre-ondelettes \\f nA dudit espace W^, 
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dans laquelle on exprime lesdits coefficients d' image. 

En d'autres termes, on remplace l'etape d' envoi des coefficients du 
maillage raffine par la prise en compte des coefficients du maillage non raffine 
(niveau precedent) auquel on rajoute 1' envoi des coefficients attribues aux 
5 nouveaux noeuds du maillage raffine. Le traitement de ces dernier s permet 
egalement d'optimiser la position des coefficients des niveaux precedents, de 
fagon a Stre optimal a chaque niveau. 

Un avantage d'une telle methode tient dans la meilleure repartition 
statistique de la valeur des coefficients permettant une diminution du cout de 
10 codage. Avec une simple transmission des coefficients en partant des sommets les 
moins raffines la nouvelle methode permet d'effectuer une premiere 
representation de l'image a partir de ces seuls coefficients optimale pour ce niveau 
de resolution. Dans les methodes connues en revanche, les premiers coefficients 
correspondent a la valeur pour ces sommets optimisee pour le niveau de 
15 raffinement le plus fin du maillage (on a done sous-optimalite). 

De fagon preferentielle, lesdites fonctions de pre-ondelettes sont des 
fonctions d'ondelettes, orthogonales entre elles. 

Avantageusement, chacune desdites fonctions de pre-ondelettes possede 
un support etroit, limite a un nombre predefini de sommets dudit maillage situe au 
20 voisinage d'un sommet de reference pour ladite fonction. 

Selon un premier mode de mise en oeuvre de Finvention, lesdites fonctions 
de pre-ondelettes sont monodimensionnelles. Dans ce cas, chacune desdites 
fonctions de pre-ondelettes est avantageusement determinee en tenant compte de 
la position dans ledit maillage d'au moins une arete portant un nouveau sommet 
25 dudit maillage, auquel ladite fonction sera affectee. 

Ainsi, la base cpi (p+1) peut notamment etre obtenue a Faide de Tequation 
suivante : 

30 ou Vi (p) est Tensemble des sommets voisins du sommet j dans le maillage de 
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niveau p. 

Selon un autre mode de realisation de l'invention, lesdites fonctions de 
pre-ondelettes sont multidimensionnelles. II peut notamment s'agir de fonctions 
de type "box splines". 

5 De fa§on avantageuse, la determination des coefficients d'image a un 

niveau de maillage n repose sur la resolution d'un systeme lineaire : 

( y(p-D^ 

avec : ^ 

v / 

ou X est le vecteur de coordonnees d'approximation de ladite image source cf p) 9 
10 W^ -0 est la matrice de passage de la base (cp^) de l'espace V p _! a la base de 
fonctions de pre-ondelettes de l'espace W p _! et V^" 1 * est la matrice de l'espace V^. 

Cette approche s'avere tres efficace, en particulier lorsque la resolution 
dudit systeme lineaire met en ceuvre un algorithme du type « matrices creuses » 
ou du type "matrices profils". 
15 Selon un aspect avantageux de l'invention, le procede de codage de 

F invention comprend : 

une etape prealable de determination de bases de fonction cp n et \|/ n 
de reference, pour une structure d'image predeterminee ; 
une etape de codage damages, utilisant systematiquement lesdites 
20 bases de fonction cp n et \|f n de reference, pour toute image ay ant 

ladite structure d'image predeterminee. 
Le codage d'un ensemble d'images de meme type est ainsi 
particulierement simplifie. 

Lors de ladite etape prealable, on determine avantageusement au moins 
25 deux jeux de bases de fonction <p n et \|/ n de reference, correspondant chacune a une 
structure d'image predeterminee particuliere. 

L'invention concerne egalement un procede de construction de bases de 
fonctions pour le codage hierarchique d'images source mettant en oeuvre un 
maillage hierarchique definissant au moins deux espaces emboites, correspondant 
30 chacun a un niveau n de decomposition dudit maillage. Selon l'invention, pour au 
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moins une structure d'image predeterminee, on associe a au moins niveau n de 
decomposition, et preferentiellement a tous les niveaux (a F exception du premier 
niveau de decomposition), un espace mathematique de representation d'une image 
V n tel que V n = 0 W n . l9 ou : 
5 - V n .! est 1' espace associe au niveau n-1 de decomposition ; 

- W n _! est un espace orthogonal a V n . l9 

et on associe audit espace V n d'une part une base de fonctions affines par 
morceaux cp n , et d' autre part une base formee de la combinaison de deux bases 
orthogonales : 

10 - une base de fonctions affines par morceaux <p nA dudit espace V n4 ; 

- une base de fonctions de pre-ondelettes \|/ n-1 dudit espace W n . l9 

dans laquelle on peut exprimer des coefficients d'image pour toute image ayant 
ladite structure d r image, permettant de produire une image reconstruite, 
representative de ladite image source, avec une qualite de restitution optimisee 
15 pour ledit niveau n de decomposition. 

L' invention concerne egalement les dispositifs de codage damages, 
mettant en oeuvre le procede de codage et/ou de construction de bases decrits ci- 
dessus. 

Ainsi, F invention concerne un dispositif de codage d'au moins une image 
20 source mettant en oeuvre un maillage hi&rarchique definissant au moins deux 
espaces emboites, correspondant chacun a un niveau n de decomposition dudit 
maillage, dans lequel a au moins un niveau n de decomposition (a l'exception du 
premier niveau de decomposition), il associe uniquement des coefficients d'image 
exprimes dans une base de fonctions definie dans un espace orthogonal a Fespace 
25 correspondant au niveau precedent n-1 de decomposition, lesdites fonctions etant 
choisies de fa§on que lesdits coefficients d'image permettent d'optimiser pour 
ledit niveau n de decomposition les informations deja transmises pour le niveau 
precedent n-1 de decomposition, de fagon a produire une image reconstruite, 
representative de ladite image source, avec une qualite de restitution optimisee 
30 pour ledit niveau n de decomposition. 
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L' invention concerne encore un precede de decodage d'images codees a 
l'aide du precede de codage et/ou de construction de bases decrits ci-dessus. 

Un tel precede de decodage decode des images codees sous la forme de 
coefficients d'image obtenus par la mise en oeuvre d'un maillage hierarchique 
5 definissant au moins deux espaces emboites, correspondant chacun k un niveau n 
de decomposition dudit maillage, selon un codage associant, a au moins un niveau 
n de decomposition (a F exception du premier niveau de decomposition), 
uniquement des coefficients d'image exprimes dans une base de fonctions definie 
dans un espace orthogonal a Fespace correspondant au niveau precedent n-1 de 
10 decomposition, lesdites fonctions etant choisies de fa?on que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
informations deja transmises pour le niveau precedent n-1 de decomposition, de 
fagon a produire une image reconstruite, representative de ladite image source, 
avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 
15 Avantageusement, un tel precede de decodage comprend : 

une etape prealable de reception et/ou de stockage de bases de 
fonctions de reference, pour au moins une structure d'image 
predeterminee ; 

une etape de decodage d'images, utilisant les bases de fonctions de 
20 reference correspondant a la structure d'image de l'image a 

decoder. 

A nouveau, cela permet de simplifier les traitements et limiter la quantite 
de donnees a transmettre, pour une serie d'images. 

L' invention concerne encore un dispositif de decodage d'images codees, 
25 caracterise en ce qu'il met en oeuvre le precede de decodage decrit ci-dessus. 

L' invention concerne egalement un programme informatique de codage 
et/ou de decodage d'images, caracterise en ce qu'il comprend des instructions de 
programme permettant la mise en oeuvre du precede de codage et/ou du precede 
de decodage decrits ci-dessus. 
30 L'invention concerne encore un signal de donnees d'images, comprend une 
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premiere partie comprenant des donnees representatives d'au moins un jeu de 
bases de fonctions de reference, pour au moins une structure d'image 
predeterminee, et une deuxidme partie comprenant des donnees representatives 
d'au moins deux images codees a l'aide d'un desdits jeux de bases de fonctions de 
5 reference. 

Avantageusement, ladite premiere partie comprend au moins deux jeux de 
bases de fonctions de reference, correspondant a des structures d'images 
distinctes, et en ce que chaque image de ladite deuxieme partie comprend une 
information permettant la selection de Tun desdits jeux de bases de fonctions de 
10 reference. 

De fa§on preferentielle, les donnees de ladite premiere partie et/ou de 
ladite deuxieme partie sont organisees de fa?on a permettre une reconstruction 
d'image progressive et/ou avec un niveau de qualite de reconstruction 
predetermine. 

15 Selon une caracteristique preferentielle, lesdites images sont codees selon 

un precede mettant en oeuvre un maillage hierarchique definissant au moins deux 
espaces emboites, correspondant chacun a un niveau n de decomposition dudit 
maillage, lesdites fonctions de reference etant choisies de fa§on que lesdits 
coefficients d'image permettent d'optimiser pour ledit niveau n de decomposition 

20 les informations deja transmises pour le niveau precedent n-1 de decomposition, 
le codage des images associant, a tout niveau n de decomposition (a 1' exception 
du premier niveau de decomposition), uniquement des coefficients d'image 
exprimes dans une base de fonctions definie dans un espace orthogonal a Fespace 
correspondant au niveau precedent n-1 de decomposition, de fagon a produire une 

25 image reconstruite, representative d'une image source, avec une qualite de 
restitution optimisee pour ledit niveau n de decomposition. 

L' invention trouve des applications dans de nombreux domaines, et en particulier 
pour le codage d'images fixes et le codage damages intra d'un flux d'images 
animees. 

30 D'autres caracteristiques et avantages de F invention apparaitront a la 
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lecture de la description suivante de modes de realisation preferentiels de 
Finvention, donnes a titre de simples exemples illustratifs et non limitatifs, et des 
dessins annexes, parmi lesquels : 

- la figure 1 illustre la fagon dont on determine les coordonnees d'un 
5 point P ; 

- les figures 2 A a 2D presentent les differentes positions d'une arete 
dans un maillage, selon l'approche de Quark, discute en annexe 1 ; 

la figure 3 est un schema synoptique simplifie du principe general d'un 
mode de realisation de 1' invention ; 
10 - la figure 4 illustre un k-disque pris en compte pour le calcul des 

fonctions de pre-ondelettes, dans un exemple de mise en oeuvre de 
F invention ; 

la figure 5 presente un exemple de structure de signal exploitant le 

principe illustre en figure 3. 
15 L'invention repose done notamment sur Futilisation d' ondelettes, ou de 

pre-ondelettes, specifiques, basee sur Forthogonalisation des bases 
complementaires. 

On peut envisager de maniere sirnilaire de Fappliquer aux ondelettes de 
type box-splines [2], Les premieres sont des preondelettes sur maillage 
20 triangulaire dont les coefficients significatifs sont peu nombreux. Les ondelettes 
box-splines sont egalement construites a partir de maillages et represented un 
moyen de construire des ondelettes bidimensionnelles qui ne soient pas basees sur 
le produit tensoriel de bases d'ondelettes monodimensionnelles. 

La methode proposee constitue une amelioration de la methode precedente 
25 de codage par maillages developpee dans [7], Elle concerne le codage d' images 
fixes mais peut aussi s'appliquer au codage des images intra dans le cadre de 
compression video en particulier les codages utilisant les maillages conjointement 
pour F approximation de surface intra et F estimation de mouvement inter. 

L'approche peut done s'appliquer de maniere generale a la transmission et 
30 au stockage d'images. 
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Parmi les applications dans lesquelles le precede peut etre utilise, on pent 
done citer : 

les applications liees a la teleconference ; 

la transmission progressive d' images naturelles sur un reseau 
5 distribue ; 

codage des images intra dans un schema de compression video base 
sur les maillages et application aux VOP. 
La figure 3 illustre done, de fagon simplifiee, le principe general de 
l'invention. 

10 On considere un maillage emboite (311), de type connu en soi, definissant 

une pluralite d'espaces V n , correspondant chacun a un niveau du maillage. 

Pour chaque de ces espaces, on recherche des bases orthogonales (312), 
selon la technique decrite en detail par la suite. On en deduit (313) des matrices de 
codage, permettant de determiner des coefficients d' image optimises pour chacun 

15 des niveaux du maillage. 

II est a noter que ces operations sont ind6pendantes du contenu de 1' image. 
Elles sont simplement liees au maillage utilise. II peut done avantageusement 
s'agir d'une etape preliminaire 31, produisant des matrices de codage et de 
decodage qui pourront ensuite etre utilisees pour une pluralite d' images, des lors 

20 qu' elles sont associees au meme maillage. 

On prevoit done une etape 32 de codage proprement dit, qui peut etre 
repetee plusieurs fois (321) sans revenir sur le calcul des bases orthogonales. On 
dispose ainsi d'un procede relativement simple, les matrices ay ant ete calculees 
une fois pour toutes. 

25 Pour chaque image source re§ue (322), qui peut Stre notamment une image 

fixe ou une image intra d'un flux video, on calcule (323) les coefficients d'image 
correspondant a chaque niveau du maillage, a l'aide des matrices, puis on 
transmet (324), ou on stocke, les coefficients obtenus. 

Le decodage 33 effectue les operations inverses. II dispose des matrices de 

30 decodage, qui peuvent etre implantees definitivement en memoire ou qui sont 
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regues (331), par exemple au debut d'un flux video, pour initialiser le decodeur. 

A reception (332) des coefficients codes, le decodeur reconstruit 
progressivement (333) chaque image jusqu'au niveau de qualite souhaitee. 

Le fait que les memes bases peuvent etre utilisees pour plusieurs images 
5 permet de definir une structure avantageuse de signal, telle qu'illustr6e en figure 5. 

Le signal comprend une premiere partie 51 ^initialisation, qui contient les 
bases de fonctions determinees une fois pour toutes, et une seconde partie 52 de 
donnees d'images, comprenant des coefficients determines a l'aide desdites bases. 
On obtient ainsi un debit reduit, la premiere partie 51 etant transmise une seule 
10 fois, en debut de sequence. 

Lorsqu f il y a plusieurs jeux de bases de fonctions 511, 
512,...correspondant a autant de types damages, on prevoit un systeme 
d'identification de celles-ci. Pour chaque image ou serie d'images, ou simplement 
lorsqu'un changement est necessaire, on transmet en preambule des donnees 521 
15 d'une image un identifiant 521 1 indiquant le jeu de bases de fonctions a utiliser. 

On peut prevoir des flux de transmission distincts pour les parties 51 et 52. 
On peut egalement prevoir au moins deux flux pour la partie 52 : un premier flux 
correspondant a un ou plusieurs niveaux "grossiers", et plusieurs autres flux 
correspondant a des niveaux plus raffines. Selon les besoins, le debit et/ou les 
20 moyens disponibles, on choisit les flux utilises. 

On peut egalement prevoir que la premiere partie 51 est connue une fois 
pour toutes des terminaux (par exemple sous la forme d'un CD-ROM). Dans ce 
cas, seule la partie 52 est transmise, ce qui correspond a un debit tres reduit. 

Le meme principe s'applique egalement, bien star, aux supports de donnees 
25 stockant des images ainsi codecs. 

On decrit maintenant plus en detail la mise en ceuvre d'un mode de 
realisation de Tinvention. 

/. Definition des nouveaux coefficients 

On considere ici que Ton dispose d'espaces emboites V 0 ,V lt . . .,V n , associes 
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a des maillages triangulares egalement emboftes. La dimension de ces espaces 
correspond au nombre de sommets du maillage considere. On suppose ici que Ton 
dispose, pour chaque espace V p d'une base cp/ p) ou chaque fonction est associee au 
sommet i de maniere univoque. En general on utilise des bases de fonctions dont 
5 la restriction a chaque triangle est polynomiale. 

Dans l'exemple decrit par la suite, par souci de simplification, on considere 
qu'un espace V n est associe a chacun des niveaux de decomposition du maillage. 
L'invention s'applique bien sur egalement au cas ou on associe des espaces 
emboites V m a V t aux seuls niveaux de decomposition consecutifs m a I du 
10 maillage, ou m et / sont deux entiers quelconques. Par exemple, le principe de 
Tinvention decrit dans la suite de ce document peut n'etre mis en oeuvre que pour 
les niveaux de decomposition 3 a 12 du maillage, associes a des espaces emboites 
V 3 aV 12 . 

La base affine presentee en annexe 2 est un cas particulier de telles bases 
15 ou les polynomes consideres sont de degre 1 c'est a dire que les espaces vectoriels 
dans lesquels on se place sont limites aux fonctions affines par morceaux et 
continues globalement. De maniere a utiliser les proprietes de la multiresolution, 
on cherche des bases qui soient reliees d'un niveau a un autre par une equation de 
raff inement de type : 

20 

(6)^= 2>^> 

On dispose grace a la hierarchie de maillages d'une suite d' espaces 
emboites, chacun correspondant a une representation et a une qualite de Timage 
donndes. Si Ton dispose des espaces V n et V n+1 on s'interesse ici a trouver une 
25 base de F espace W n tel que V n ®W n = V n+1 ou le ® designe la somme orthogonale 
des 2 espaces vectoriels. On sait que Ton peut done avoir l'ecriture suivante pour 
P approximation de Pimage I : 



WO 02/100111 



14 



PCT/FR02/01937 



/=i P =i i=\ 

ou l'indice p decrit P ensemble des differentes resolutions et ou N p designe la 
dimension de W p . 

Les fonctions \\f ip \ sont choisies de maniere a constituer une base de 
5 l'espace W p . Par ailleurs, dim W p = dim V p+1 - dim V p , done le nombre de 
fonctions \|/ (p \ est egal au nombre de nouveaux sommets pour la resolution p. On 
peut done associer de maniere bijective les fonctions de base de W p avec les 
nouveaux sommets. 

On appelle ces fonctions preondelettes. Dans le cas particulier ou les 
10 fonctions de base d'un meme niveau sont orthogonales entre elles on les appellera 
ondelettes. 

De maniere pratique, les obstacles rencontres dans P application de la 
methode proposee proviennent de la difficulty d'exhiber les bases (\|/ (p) i) dont les 
fonctions possedent des supports limites a un nombre fini et le plus faible possible 

15 de triangles contigus au sommet de reference. Ainsi on ne peut se contenter 
d'orthogonaliser la base initiale de maniere arbitraire. On donne dans les sections 
suivantes des exemples pratiques d' application ou la methode peut s'appliquer. 

Une fois que Ton dispose d'une telle decomposition des approximations de 
F image, on peut coder les coefficients obtenus selon les m^thodes classiques. 

20 Notons que dans ce cas l'orthogonalite entre deux niveaux de resolution 
successifs invalide un codage differentiel : les coefficients obtenus sont 
directement quantifies et codes par le biais d'un codeur arithmetique. 

2. Calcul pratique des nouveaux coefficients 

Pour simplifies le calcul est detaille dans le cas ou nous disposons 
25 simplement de deux espaces successifs V 0 et On cherche done une base de 
l'espace W 0 . Comme chaque fonction 0|/ (0) i) de cette base appartient a l'espace Vj 
et que V 0 ®W 0 = V x on sait qu'il existe des coefficients w a tels que : 
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(1) 



10 



15 



La premiere etape de ralgorithme consiste done a determiner la matrice 



En effet cette matrice constitue la matrice de passage de la base (<pi (1) ) de 
l'espace V x a la base de preondelettes de l'espace W 0 . On remarque que ces 
calculs sont independants du contenu de l'image, et ne dependent en fait que de la 
forme du maillage et du type de base choisie. II est done possible d'effectuer un 
precalcul pour les maillages classiquement utilises. Ce precalcul depend de la 
maniere dont on construit les bases de preondelettes. La sous-section 5.4 et la 
section 6 montrent des exemples pratiques de telles constructions. De meme on 
note W(p) la meme matrice pour les niveaux p et p+1. 

Decrivons maintenant le calcul effectif des coefficients a partir des 
matrices de passage precedentes. Si pour une image donnee on dispose de p 
niveaux de resolution : 



i=\ 

On peut calculer pour chaque niveau les coefficients (P(p)i) qui 
correspondent aux coefficients de 1' equation (7) en resolvant un systeme lineaire. 
En effet on a : 



W< 0) =(wA 




(10) v ip ~ l) a ip " l) =a ip) 



(11) W^P^'V =a ip) 



II suffit done pour obtenir les coefficients voulus de resoudre le systeme 
lineaire 
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(12) A {p ~ l) X = a {p) 



avec 



Le vecteur X etant compose des coordonnees d' approximation de 1' image 
5 dans la nouvelle base. 

Parmi les methodes de resolution effective de tels systemes lineaires, on 
peut distinguer differentes methodes parmi lesquelles les methodes iteratives. On 
s'est plutot oriente vers des methodes de resolution basees sur 1' utilisation des 
proprietes des matrices creuses. Naturellement cette methode n'est applicable 
10 dans la pratique que dans la mesure ou les matrices correspondantes sont 
suffisamment creuses pour permettre une complexity en temps de calcul 
acceptable grace a 1' utilisation des matrices profils. 

Dans notre cas precis, cette propriete correspond a la limitation du support 
des fonctions preondelettes utilisees a un nombre suffisamment faible de 
15 sommets. Ces methodes de resolution de systemes lineaires constituent une 
premiere possibility. II peut cependant exister d'autres solutions susceptibles 
d'accelerer les temps de calcul et evitant la resolution explicite de Tintegralite du 
systeme lineaire. 

Nous allons illustrer par quelques exemples comment il est possible 
20 d'exhiber des bases de preondelettes sur des maillages. 

3 . 1 determination des bases 

On donne ici Fexemple de preondelettes basees sur les fonctions elements 
finis affines de Lagrange. On definit le k-disque D ik (p) autour d'un sommet i du 
maillage comme T ensemble des sommets que Ton peut relier par un chemin 
25 constitue de moins de k aretes du maillage du niveau p (voir figure 3 ). On 
cherche ici simplement a construire la base de preondelettes comme etant la 
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projection orthogonale de (<p (p+1) j) sur W p a support limite a un k-disque autour du 
sommet i, ou i decrit r ensemble des noeuds du niveau p+1 qui n'existaient pas au 
niveau p. On exprime les fonctions (\)/ (p) i) de la maniere suivantepour les i 
correspondant aux sommets qui appartiennent au niveau p mais n'existaient pas 
5 au niveau p-1: 



00 



(13)?^=^ +1) - 2^} 

Pour determiner les coefficients de la matrice Q=qij on doit ecrire les 
conditions d'orthogonalite de \\f (p \ avec V p . On ecrit done : 

\/k = l,...,N p ,(V?\<pi P) ) = 0 

10 soit encore (pour i nouveau sommet et k ancien sommet): 



Grace a 1' equation de raffinement, on connait 1' expression des fonctions de 
V p en fonction de fonctions de V p+1 , On en deduit que Q est solution d'un systeme 
15 lineaire qui ne depend que des produits scalaires de fonctions de meme niveau. On 
sait calculer ces produits scalaires en fonction de la mesure des triangles [6]. 



2 jeV .tP) 



On est maintenant en etat de deduire les coefficients de la matrice W en se 
20 servant de l'equation de raffinement (4). Dans la sous-section suivante, on 
explicite les procedures utiles pour le calcul et le codage de Fimage. 
3.2 Calcul et codage 
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Supposons que Ton veuille coder une image sur laquelle on dispose de 
deux niveaux de maillage. La premiere etape de ralgorithme consiste a effectuer 
par la methode decrite dans [7] une optimisation des valeurs nodales sur la base 
d'elements finis affines cp (p) i9 ou i decrit l'ensemble des sommets du niveau 1. On 
5 obtient done la decomposition donnee par l'equation (3). 

On rappelle que la minimisation de l'erreur quadratique definie dans 
l'equation (2) revient a la resolution d'un systeme lineaire que Ton resout par des 
methodes profils [7]. On ecrit ce systeme lineaire de la maniere suivante : 

E= 2 (Kx,y)-^ai P Vn P) (^y)f 



10 On consid&re done l'image definie par ses valeurs pixeliques I(x,y). On 

cherche a minimiser E, ou Q represente l'ensemble des pixels de l'image et Vp 
l'ensemble des sommets du niveau de resolution p. On obtient le minimum en 
ecrivant les Np equations : 

dE A 1 AT 



dai p) 



15 Pour un sommet n donne on note S(n) le support de la fonction (p(p)n, on a 

alors les Np equations suivantes : 

S E Qn & y)<Pk (*> yy*F* = S J ( x ' yypn y) 

(x,y)eS(n) keV^ (x f y)ES(n) 



On doit done resoudre done un systeme linaire Ax=B ou Ton calcule la 
matrice A et le vecteur B de la maniere suivante, apres avoir initialise A et B a 0 : 
20 Pour i=1,...,N T (nombre de triangles pour le maillage considere) 
Pour (x,y) dans T, (parcours des pixels du triangle T,) 
Pour j=1 a 3 (sommets de T,) 
Calcul de f, = cp mG) (x,y) ou m(j) est I'indice du j omo sommet de T, 
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Pour (k=1 a 3, j =1 a 3) (rem. Si (k,j) est fait on ne fait pas Q,k)) 

A m ( k ) |m (j) = A m(k)jm{J) + f k * f j 
Pour (|= 1 a 3) 
B ma) = B m(J) + fj * l(x,y) 
5 Fin pour j 

Fin pour i 

Algorithme 1 : Optimisation des valeurs nodales 

La matrice A etant symetrique definie positive et presentant un nombre 
important de valeurs nulles, on resout alors le systeme lineaire par une methode 
10 profit. 

Maintenons que Ton dispose des coefficients a determinons les 
coefficients sur la base des preondelettes. De l'equation (15) on deduit un systeme 
lineaire permettant de decomposer 1' image sous la forme (7) pour P=p. Decrivons 
le passage entre deux niveaux c'est-a-dire le passage de l'ecriture : 

15 

aTecritxire : 

z=l k=p i=l 

A chaque etape on cherche done a resoudre un systeme lineaire Cx=oc ou 
les coefficients de la matrice C sont les coefficients de synthese. On initialise la 
20 matrice C a 0. Dans le cas present, on va determiner la matrice C en utilisant les 
equations (4) et (13). Dans un premier temps, donnons 1'algorithme qui permet de 
calculer les coefficients qij de l'equation (13). Pour cela on deduit un systeme des 
equations (14) Dq=E. 

Pour simplifier on considere ici le k-disque egal a tout le maillage. On 
25 initialise les matrices D et E a 0. Avec : 
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E lk =«prM p) > 

B u =«pf p) M P+1) > 

Si Ton choisit un maillage pour lequel tous les triangles d'un meme niveau 
ont la meme aire :on a les relations suivantes qui permettent un calcul aise des 
5 coefficients precedents [6] : 
Pour i^j : 

<<p?\qffb>=l YAiT) 

Pour i = j : 

TeS(i)nS(j) 

Ce nouveau systeme lineaire est encore un systeme a matrice suffisamment 
creuse ce qui nous autorise a utiliser une methode adaptee. Cependant la matrice 
n'etant pas symetrique positive, Pour appliquer les methodes profil, on remplace 
done le systeme lineaire Cx-a par le systeme equivalent tCCx = tCa. On obtient 
ainsi un systeme lineaire dont la matrice correspondante est bien symetrique 
definie positive (en effet C est inversible). Ce systeme est moins creux, mais la 
matrice correspondante possede encore un profil acceptable. 

Lorsque Ton dispose de plus de deux niveaux de maillages, on itere le 
processus precedemment decrit. Une fois que Ton dispose de cette decomposition 
sur 1' ensemble des niveaux de resolution, e'est-a-dire des coefficients a(0)i et 
p(p)i de P equation (7), on effectue une quantification scalaire et on transmet les 
valeurs quantifiees correspondant a un sommet selectionne dans Parbre de 
triangles raffines. 



10 



15 



20 
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La phase de decodage s'effectue en synth6tisant les coefficients, cette fois- 
ci directement avec la matrice C et non son inverse. On obtient de cette maniere 
des coefficients qui permettent d'obtenir une reconstruction globale de l'image 
par simple interpolation. 
5 4. deuxieme exemple 

On donne ici un autre exemple de schema de codage d' images fixes 
comme application de la presente methode. On utilise ici les preondelettes basees 
aretes decrites dans Tannexe 1. La phase d' optimisation des valeurs nodales est 
commune avec le premier exemple. 
10 Dans ce cas le calcul des coefficients d'ondelettes s'effectue egalement a 

l'aide de la resolution d'un systeme lineaire Cx=a. Pour determiner les 
coefficients de la matrice C on applique alors 1' algorithms suivant : 
Pour i=1,...,N p 
C„ = 1 

15 Pour j voisin de i dans le maillage 

C u = 1 
Fin pourj 
Fin pour i 

Pour i= Np+1 N p+1 

20 C„ = U determine par la position de I'arete i correspondante 

dj ^ O est calcule pour les cas determines dans sous-section 
precedente 
Fin pour i 

Algorithme 2 : Remplissage de la matrice de passage C. 
25 Notons que cette methode ne pose pas d' implementation sur les bords. En 

effet le mode de construction des preondelettes a tenu compte des differentes 
configurations possibles pour un maillage borne quelconque. 
5. application avec des "box splines " 
5A De scription de la methode 
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On peut egalement appliquer la methode sur d'autres types de construction 
d'ondelettes ou de preondelettes. On peut ainsi construire des bases explicites 
d'ondelettes a partir d'espaces emboites sur lesquels on a construit une multi- 
resolution. Une methode generate de construction de bases d'ondelettes multi- 
5 dimensionnelles est developpee dans [3]. L'interet d'utiliser de telles ondelettes 
par rapport aux ondelettes construites par produit tensoriel reside dans la 
meilleure adaptation des fonctions de base au cadre multi-dimensionnel. 

Dans [2] ce principe est repris et on donne de plus une methode permettant 
de determiner une base de preondelettes a support compact qui seront done 
10 directement utilisables dans la pratique. On se propose ici d' appliquer cette 
decomposition en ondelettes box-splines 

5.2 Algorithme 

Des exemples de fonctions box-splines sont donnees dans [2]. Les 
coefficients explicites des fonctions de la base y sont donnes. On peut done, 
15 comme precedemment en deduire une matrice de passage W qui nous permet de 
determiner les coefficients ondelettes correspondant. Les algorithmes s'ecrivent 
exactement de la fagon suivante : par exemple les coefficients peuvent etre deduits 
directement des exemples explicites dans [2], 
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ANNEXE 1 

Une arete est dite interieure si elle dflimite deux triangles du maillage, 
sinon elle est dite exterieure. L'ardte peut Stre interieure avec ses deux sommets 
interieurs(a), interieure avec Tun des deux sommets exterieur (b), interieure avee 
5 ses deux sommets exterieurs (c) ou encore elle meme exterieure (d). Les differents 
cas sont montres figure2. Une telle base de representation presente done 
Favantage d'etre adaptee a la multiresolution et d'gtre utilisable facilement gr&ce 
a la petite taille du support utilise. 

Rappelons 1' expression explicite des preondelettes de Quak dans les 
10 differents cas. Considerons le premier cas. Soit u le nouveau sommet considere, 
soient a x et a 2 les deux sommets et s et t leurs connectivites respectives au niveau 
de maillage 1, on a alors une preondelette possible \|/ u (0) en ecrivant : 

(16) ¥ ™ = Ucpf + A x q>™ + A 2 (p% + + 

ou bj decrit les sommets connexes a a 2 et differents de u et Ci les sommets 
15 connexes a a 2 et differents de u. Pour un sommet a donne la fonction (p a designe ici 
F element fini affine associe au sommet a pour le maillage correspondant. 
U expression des coefficients est donnge par : 

4 = -^ 
3V2T 



it 

U = -—(A 1 +A 2 ) + 



1 A M+Af-' 

B. = A + ; 

14 1 \-X s 

r 1 . 



14 1 \-X 
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avec : 



2 2 

Dans le deuxieme cas, on suppose que l'une des deux extremites de 1' arete 
est un sommet exterieur du maillage. On note alors s x le nombre de sommets 
5 voisins de a 2 et situes a gauche de 1' arete consideree et s 2 le nombre de sommets 
voisins de a x et situes a droite de 1' arete. Les coefficients correspondant sont 
respectivement notes B^et B 2 ^ On a alors l'expression suivante des coefficients : 

3->/2T 



4=" 

A 2 = 



2(s,+5 2 ) 

3V2T 



2t 

B) =-—4+ XS2+ \ S \ (A'"~ J| +A^'-'') 
14 A - * 1- * 2 -A s ' +S2 

/ 1/1 



14 * 

Dans le troisieme cas on utilise les memes notations en distinguant 
10 egalement les t x sommets voisins de b x et situes a gauche de F arete consideree et ^ 
le nombre de sommets voisins de b x et situes a droite de F arete. 
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3V21 

A = 



A 



3V21 



2(r 1 +/ 2 ) 

14 1 2 * AT*'* 2 -% S,+S2 A - ' 1- ' 2 -A' 1+ ' 2 

La > ^ + ^~ S2 

14 1 + A~*~* 2 -A* 1+J2 



B) =-—4 + - . , s+s (A^ 1 +A S '"'') 



1 3 s i -1- a — ^"t 

B ' 77 Ai + _ I*** ^ +A } 



14 1 A'^-A^ 

c; =-—4 + .f, +A ~l (A' Wl +A"-'') 

14 2 A 1 2 - A 1 2 
J_ 
14 



1 , A''+A-'' 



Enfin dans le dernier cas, ou 1' arete est exterieure on a : 

1 25 

2 2/ 

C/— 1 (A + A )+ 2(A- + AQ 2(A- + A) 



5 
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On dispose done dans tous les cas de figures de P expression de la 
preondelette correspondant a un nouveau sommet en fonction de la base des 
fonctions affines du niveau le plus fin. 
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ANNEXE 2 



1. Approximation de surfaces 

Pour simplifier, on considere une image en niveaux de gris (luminance). 
Pour les bandes de chrominance, on peut appliquer une methode similaire. On 
5 peut done considerer F image comme 6tant une representation discretisee d'une 
surface parametrique. On s'interesse ici a la representation de cette surface par un 
maillage. 

On suppose que Ton dispose d'un maillage triangulaire quelconque 
recouvrant 1' ensemble du domaine de T image. On utilise alors les fonctions de 
10 base C° de Lagrange centrees sur un sommet. On obtient ainsi une approximation 
/ de V image I sous la forme : 



ou N designe le nombre de sommets du maillage, ^ est la fonction affine prenant 
15 la valeur 1 sur le sommet i et telle que la valeur en un point d'un des triangles 
dont i est le sommet soit egale a la coordonnee barycentrique de ce point par 
rapport au sommet i (voit figure 1). Cette fonction est done nulle en dehors des 
faces triangulares dont Tun des sommets est le sommet i. Les coefficients oq sont 
calcules de maniere a minimiser l'erreur quadratique: 



N 



(l)/ = ^a.<p /9 V/,a.ei? 




20 



sur le support de 1' image. 
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2. Hierarchie de maillages 

La methode utilise une hierarchie de maillages triangulares emboftes. On 
dispose initialement d'un maillage triangulaire. On se dote ensuite d'une regie de 
subdivision des maillages. On peut par exemple obtenir les maillages suivants en 
5 inserant au milieu de chaque arete du maillage un nouveau sommet. 

Chaque triangle est dans ce cas divise en 4 nouveaux triangles. On utilise 
cette representation de F image de manure a disposer d' approximations 
successives de F image, chacune correspondant a une resolution donnee. On 
exploite ensuite cette hierarchie en subdivisant uniquement les triangles dont 
10 Ferreur quadratique de reconstruction est superieure a un certain seuil. Le valeur 
de ce seuil determine en partie la qualite de reconstruction desiree. On obtient 
ainsi un arbre de subdivision de la hierarchie qui doit etre transmis au codeur et 
permet au decodeur de determiner quels sont les triangles raf fines. 

3. Codage des coefficients 

15 On transmet dans tous les cas des valeurs pour le niveau de base (le plus 

grossier de la hierarchie). En effet ces valeurs sont necessaires pour une 
reconstruction minimale de Fimage. Au niveau 0 on a : 



Ou les N 0 premiers sommets correspondent aux sommets du maillage de niveau 1 
qui existaient deja dans le maillage de niveau 0 et les N x N 0 sommets du deuxieme 




20 Au niveau 1 on peut done ecrire : 
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terme correspondent aux sommets du niveau 1 qui n'existaient pas au niveau 0. A 
un niveau p quelconque on dispose done de la representation : 

(3) / = fX'V,^ = %a\^ + %j p ty p) 

5 Le premier terme de la somme correspond aux N p sommets du niveau de 

resolution p et le deuxieme terme de la somme correspond aux nouveaux sommets 
inseres au niveau p+1. Pour les N p premiers sommets on a : 

10 Ou Vj*' represente 1' ensemble des sommets voisins du sommet i dans le maillage 
P- 

On peut done exprimer la representation de 1 'image sous la forme 
diff6rentielle suivante : 

;=1 p=0 i=l 

15 4. Remarque : 

De maniere a homogeneiser les coefficients ainsi obtenus, on retranche aux 
coefficients du premier niveau leur moyenne. Cette representation permet 
d' exploiter les redondances spatiales liees a la correlation statistique entre pixels 
proches dans l'image. La compression est effectuee par le biais d'une 

20 quantification de ces coefficients suivie d'un codage arithm6tique adaptatif. Par 
ailleurs les coefficients n'appartenant a aucun des triangles selectionnes dans 
Tetape de raffinement adaptatif ne sont pas codes ; en effet ne pas raffiner un 
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triangle revient a considerer que les coefficients des sommets correspondant aux 
resolutions plus fines a l'interieur de ce triangle sont les valeurs interpolees, ce qui 
est equivalent a assigner 0 au coefficient Si (p) . 

5. Interet de la methode 

5 La methode permet done d'obtenir un schema simple de codage d' images 

fixes par le biais de V utilisation d'une hierarchie de maillages reguliers emboites. 
Ceci permet d'obtenir des taux de compression efficaces. Par ailleurs un tel 
schema de codage est bien adapte a une transmission scalaire des coefficients. 
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REVENDICATIONS 

1. Procede de codage d'au moins une image source mettant en oeuvre un 
maillage hi6rarchique definissant au moins deux espaces emboftes, correspondant 
chacun a un niveau de decomposition dudit maillage, 

5 caracterise en ce qu'on associe a au moins un niveau n de decomposition (a 
l'exception du premier niveau de decomposition) un espace mathematique de 
representation d'une image V n tel que V n = V n _! ® W n _ 1? ou : 

- est Fespace associe au niveau n-1 de decomposition ; 

- W n _! est un espace orthogonal a V n „ 1? 

10 en ce qu'on associe audit espace V n d'une part une base de fonctions affines par 
morceaux <p n , et d' autre part une base formee de la combinaison de deux bases 
orthogonales : 

- une base de fonctions affines par morceaux cp^ dudit espace ; 

- une base de fonctions de pre-ondelettes \\f nA dudit espace W n . l5 

15 et en ce que, audit niveau n de decomposition (a 1' exception du premier niveau de 
decomposition), on delivre uniquement des coefficients d' image exprimes dans 
ladite base de fonctions de pre-ondelettes, 

lesdites fonctions de pre-ondelettes etant choisies de fa§on que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
20 informations deja transmises pour le niveau precedent n-1 de decomposition, 

de fagon a produire une image reconstruite, representative de ladite image source, 
avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 

2. Procede de codage selon la revendication 1, caracterise en ce que, a tous 
lesdits niveaux de decomposition (& l'exception dudit premier niveau de 

25 decomposition), on delivre uniquement des coefficients d'image exprimes dans 
ladite base de fonctions de pre-ondelettes definie dans un espace orthogonal a 
l'espace correspondant au niveau precedent n-1 de decomposition. 

3. Procede de codage selon Tune quelconque des revendications 1 et 2, 
caracterise en ce que lesdites fonctions de pre-ondelettes sont des fonctions 

30 d'ondelettes, orthogonales entre elles. 
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4. Procede de codage selon Tune quelconque des revendications 1 a 3, 
caracterise en ce que chacune desdites fonctions de pre-ondelettes possede un 
support etroit, limite a un nombre predefini de sommets dudit maillage situe au 
voisinage d'un sommet de reference pour ladite fonction. 
5 5. Procede de codage selon Tune quelconque des revendications 1 a 4, 
caracterise en ce que lesdites fonctions de pre-ondelettes sont 
monodimensionnelles. 

6. Procede de codage selon la revendication 5, chacune desdites fonctions de 
pre-ondelettes est determinee en tenant compte de la position dans ledit maillage 

10 d'au moins une arete portant un nouveau sommet dudit maillage, auquel ladite 
fonction sera affectee. 

7. Procede de codage selon Tune quelconque des revendications 1 a 6, 
caracterise en ce que la base (pi (p+1) est obtenue a Taide de Tequation suivante : 

15 <pr i)= 9^+i ^j>f 

ou Vi (p) est l'ensemble des sommets voisins du sommet j dans le maillage de 
niveau p. 

8. Procede de codage selon Tune quelconque des revendications 1 a4, 
caracterise en ce que lesdites fonctions de pre-ondelettes sont 

20 multidimensionnelles. 

9. Procede de codage selon la revendication 8, caracterise en ce que lesdites 
fonctions de pre-ondelettes sont des fonctions de type "box splines". 

10. Procede de codage selon Tune quelconque des revendications 1 a 9, 
caracterise en ce que la determination des coefficients d'image a un niveau de 

25 maillage n repose sur la resolution d'un systSme lineaire : 

A (p ~ 1) X=of p) 



avec : A^'^- 



v J 

ou X est le vecteur de coordonnees d'approximation de ladite image source oft\ 
W^" l) est la matrice de passage de la base ((pi (p l) ) de Tespace V p _ x a la base de 
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fonctions de pr6-ondelettes de 1'espace W p-1 et V*" 0 est la matrice de 1'espace V v . v 
11. Procede de codage selon la revendication 10, caracterise en ce que la 
resolution dudit systeme lineaire met en oeuvre un algorithme du type « matrices 
creuses » ou du type "matrices profils". 
5 12. Procede de codage selon Tune quelconque des revendications 1 a 11, 
caracterise en ce qu ? il comprend : 

une etape prealable de determination de bases de fonction (p n et \|/ n 

de reference, pour une structure d'image predeterminee ; 

une etape de codage d'images, utilisant systematiquement lesdites 
10 bases de fonction <p n et \|/ n de reference, pour toute image ay ant 

ladite structure d'image predeterminee. 

13. Procede de codage selon la revendication 12, caracterise en ce que, lors de 
ladite £tape prealable, on determine au moins deux jeux de bases de fonction <p n et 
\|/ n de reference, correspondant chacune a une structure d'image predeterminee 

15 particuliere. 

14. Procede de construction de bases de fonctions pour le codage hierarchique 
d'images source mettant en oeuvre un maillage hierarchique definissant au moins 
deux espaces emboites, correspondant chacun a un niveau de decomposition dudit 
maillage, 

20 caracterise en ce que, pour au moins une structure d'image predeterminee, on 
associe a au moins un niveau n de decomposition Ql 1' exception du premier niveau 
de decomposition), un espace mathematique de representation dune image V n tel 
queV^V^eW^ou: 

- V n _i est 1'espace associe au niveau n-1 de decomposition ; 
25 - W n _! est un espace orthogonal a V n „ x , 

et en ce qu'on associe audit espace V n d'une part une base de fonctions affines par 
morceaux <p n , et d' autre part une base formee de la combinaison de deux bases 
orthogonales : 

- une base de fonctions affines par morceaux cp n . x dudit espace V n _ t ; 
30 - une base de fonctions de pre-ondelettes \\f nA dudit espace W n . ls 
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dans laquelle on peut exprimer des coefficients d' image pour toute image ay ant 
ladite structure d'image, perxnettant de produire une image reconstruite, 
representative de ladite image source, avec une qualite de restitution optimisee 
pour ledit niveau n de decomposition. 
5 15. Dispositif de codage d'images, mettant en oeuvre le proc6de de codage 
et/ou de construction de bases selon Tune quelconque des revendications 1 a 14. 

16. Dispositif de codage d'au moins une image source mettant en oeuvre un 
maillage hierarchique definissant au moins deux espaces emboites, correspondant 
chacun a un niveau de decomposition dudit maillage, 

10 caracterise en ce que, a au moins un niveau n de decomposition (a F exception du 
premier niveau de decomposition), il associe un espace mathematique de 
representation d'une image V n tel que V n = V n _! ® W n . l9 ou : 

- V,,.! est F espace associe au niveau n-1 de decomposition ; 

- W n _! est un espace orthogonal a V n . l5 

15 en ce qu'il associe audit espace V n d'une part une base de fonctions af fines par 
morceaux cp n , et d' autre part une base formee de la combinaison de deux bases 
orthogonales : 

- une base de fonctions affines par morceaux (p n-1 dudit espace ; 

- une base de fonctions de pre-ondelettes \|/ n _! dudit espace W n _!, 

20 et en ce que, audit niveau n de decomposition (a F exception du premier niveau de 
decomposition), il associe uniquement des coefficients d'image exprimes dans 
ladite base de fonctions de pre-ondelettes, 

lesdites fonctions de pre-ondelettes etant choisies de fa§on que lesdits coefficients 
d' image permettent d'optimiser pour ledit niveau n de decomposition les 
25 informations deja transmises pour le niveau precedent n-1 de decomposition, 

de fa9on a produire une image reconstruite, representative de ladite image source, 
avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 

17. Dispositif de codage selon la revendication 16, caracterise en ce que, a 
tous lesdits niveaux de decomposition (a F exception dudit premier niveau de 

30 decomposition), il associe uniquement des coefficients d'image exprimes dans 
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ladite base de fonctions de pre-ondelettes. 

18. Procede de decodage d'images codees a l'aide du procede de codage et/ou 
de construction de bases selon Tune quelconque des revendications 1 a 14. 

19. Procede de decodage d'images codees sous la forme de coefficients 
5 d'image obtenus par la mise en oeuvre d ! un maillage hierarchique definissant au 

moins deux espaces emboJtes, correspondant chacun a un niveau de 
decomposition dudit maillage, 

selon un codage associant, a au moins un niveau n de decomposition (a 
T exception du premier niveau de decomposition) un espace mathematique de 
10 representation d'une image V n tel que V n = V n _j © W n _ u oix : 

- V n-1 est T espace associe au niveau n-1 de decomposition ; 

- W n _! est un espace orthogonal a V n _ 1? 

ledit codage associant egalement audit espace V n d'une part une base de fonctions 
affines par morceaux cp n , et d' autre part une base formee de la combinaison de 
15 deux bases orthogonales : 

- une base de fonctions affines par morceaux cp^ dudit espace V n .j ; 

- une base de fonctions de pre-ondelettes \\f nA dudit espace W n-1 , 

ledit codage associant en outre, audit niveau n de decomposition (a P exception du 
premier niveau de decomposition), uniquement des coefficients d'image exprimes 

20 dans ladite base de fonctions de pre-ondelettes, 

lesdites fonctions de pre-ondelettes etant choisies de fagon que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
informations deja transmises pour le niveau precedent n-1 de decomposition, 
de fagon a produire une image reconstruite, representative de ladite image source, 

25 avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 

20. Proced6 de decodage selon la revendication 19, caracterise en ce que, a 
tous lesdits niveaux de decomposition (a P exception dudit premier niveau de 
decomposition), on associe uniquement des coefficients d'image exprimes dans 
ladite base de fonctions de pre-ondelettes. 

30 21. Procede de decodage selon Pune quelconque des revendications 19 ou 20, 
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le procede de codage etant un precede selon Tune quelconque des revendications 1 

a 13, caracterise en ce qu'il comprend : 

une etape prealable de reception et/ou de stockage de bases de 
fonctions de reference, pour au moins une structure d'image 
5 predeterminee ; 

une etape de decodage d'images, utilisant les bases de fonctions de 
reference correspondant a la structure d'image de l'image a decoder. 

22. Dispositif de decodage d'images codees, caracterise en ce qu'il met en 

oeuvre le procede de Tune quelconque des revendications 19 a 21. 
10 23. Dispositif de decodage d'images codees sous la forme de coefficients 

d'image obtenus par la mise en oeuvre d'un maillage hierarchique definissant au 

moins deux espaces emboites, correspondant chacun a un niveau de 

decomposition dudit maillage, 

selon un codage associant, a au moins un niveau n de decomposition (a 
15 F exception du premier niveau de decomposition) un espace mathematique de 
representation d'une image V n tel que V n = V n-J ® W n . l9 ou : 

- V n „! est F espace associe au niveau n-1 de decomposition ; 

- W n _! est un espace orthogonal a V n . l5 

ledit codage associant egalement audit espace V n d'une part une base de fonctions 
20 affines par morceaux (p n , et d' autre part une base formee de la combinaison de 
deux bases orthogonales : 

- une base de fonctions affines par morceaux cp n . x dudit espace ; 

- une base de fonctions de pre-ondelettes \|/ n _! dudit espace W n-1 , 

ledit codage associant en outre, audit niveau n de decomposition (a F exception du 
25 premier niveau de decomposition), uniquement des coefficients d'image exprimes 
dans ladite base de fonctions de pre-ondelettes, 

lesdites fonctions de pre-ondelettes etant choisies de fa9on que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
informations deja transmises pour le niveau precedent n-1 de decomposition, 
30 de fa§on a produire une image reconstruite, representative de ladite image source, 
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avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 

24. Dispositif de decodage selon la revendication 23, carac tense en ce que, a 
tous lesdits niveaux de decomposition (k l'exception dudit premier niveau de 
decomposition), on associe uniquement des coefficients d' image exprimes dans 

5 une base de fonctions definie dans un espace orthogonal a Tespace correspondant 
au niveau precedent n-1 de decomposition. 

25. Programme informatique de codage et/ou de decodage d'images, 
caracterise en ce qu'il comprend des instructions de programme permettant la mise 
en oeuvre du procede de codage selon Tune quelconque des revendications 1 a 14 

10 et/ou du procede de decodage selon l'une quelconque des revendications 19 a 21. 

26. Signal de donnees d'images, caracterise en ce qu'il comprend une premiere 
partie comprenant des donnees representatives d'au moins un jeu de bases de 
fonctions de reference, pour au moins une structure d'image predeterminee, et une 
deuxieme partie comprenant des donnees representatives d'au moins deux images 

15 codees a l'aide d'un desdits jeux de bases de fonctions de reference. 

27. Signal de donnees d'images selon la revendication 26, caracterise en ce que 
ladite premiere partie comprend au moins deux jeux de bases de fonctions de 
reference, correspondant a des structures d'images distinctes, et en ce que chaque 
image de ladite deuxieme partie comprend une information permettant la selection 

20 de Tun desdits jeux de bases de fonctions de reference. 

28. Signal de donnees d'images selon Tune quelconque des revendications 26 
et 27, caracterise en ce que les donnees de ladite premiere partie et/ou de ladite 
deuxieme partie sont organisees de fa§on a permettre une reconstruction d'image 
progressive et/ou avec un niveau de qualite de reconstruction predetermine. 

29. Signal de donnees d'images selon Tune quelconque des revendications 26 a 
28, caracterise en ce que lesdites images sont codees selon un procede mettant en 
oeuvre un maillage hierarchique definissant au moins deux espaces emboites, 
correspondant chacun a un niveau de decomposition dudit maillage, 
lesdites fonctions de reference etant choisies de fagon que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
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informations deja transmises pour le niveau precedent n-1 de decomposition, 
le codage des images associant, a tout niveau n de decomposition (a 1' exception 
du premier niveau de decomposition) un espace mathematique de representation 
d'une image V n tel que V n = ® W n . l9 ou : 
5 - V n-1 est 1' espace associe au niveau n-1 de decomposition ; 

- est un espace orthogonal a V n _j, 

ledit codage associant egalement audit espace V n d'une part une base de fonctions 
affines par morceaux (p n , et d' autre part une base formee de la combinaison de 
deux bases orthogonales : 
10 - une base de fonctions affines par morceaux (p nA dudit espace ; 

- une base de fonctions de pre-ondelettes \y nml dudit espace W n . 1? 

ledit codage associant en outre, audit niveau n de decomposition (a l'exception du 
premier niveau de decomposition), uniquement des coefficients d' image exprimes 
dans ladite base de fonctions de pre-ondelettes, 

15 lesdites fonctions de pre-ondelettes etant choisies de fagon que lesdits coefficients 
d'image permettent d'optimiser pour ledit niveau n de decomposition les 
informations deja transmises pour le niveau precedent n-1 de decomposition, 
de fa§on a produire une image reconstruite, representative de ladite image source, 
avec une qualite de restitution optimisee pour ledit niveau n de decomposition. 

20 30. Application du procede de codage selon Tune quelconque des 
revendications 1 a 13 et/ou du procede de decodage selon Tune quelconque des 
revendications 19 a 21 au codage d'images fixes. 

31. Application du procede de codage selon Tune quelconque des 
revendications 1 a 13 et/ou du procede de decodage selon Tune quelconque des 
25 revendications 19 a 21 au codage et/ou au decodage d'images intra d'un flux 
d'images animees. 
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